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ABSTRACT 

We provide a detailed description of the Fortran code CPsuperH, a newly-developed com- 
putational package that calculates the mass spectrum and decay widths of the neutral and 
charged Higgs bosons in the Minimal Supersymmetric Standard Model with explicit CP 
violation. The program is based on recent renormalization-group-improved diagrammatic 
calculations that include dominant higher-order logarithmic and threshold corrections, b- 
quark Yukawa-coupling resummation effects and Higgs-boson pole-mass shifts. The code 
CPsuperH is self-contained (with all subroutines included), is easy and fast to run, and is 
organized to allow further theoretical developments to be easily implemented*. The fact 
that the masses and couplings of the charged and neutral Higgs bosons are computed at a 
similar high-precision level makes it an attractive tool for Tevatron, LHC and LC studies, 
also in the CP-conserving case. 



*The program may be obtained from 



http://theory.ph.man.ac.uk/~jslee/CPsuperH.html 



1 Introduction 



The quest for the still-elusive Higgs boson pDj , the missing cornerstone of the renormalizable 
Standard Model (SM), has become not only more pressing, after the completion of the LEP 
experimental programme, but also more exciting in light of the upcoming experiments 
at the upgraded Tevatron collider and the Large Hadron Collider (LHC). Indeed, direct 
searches for possible realizations of the mechanism of spontaneous electroweak symmetry 
breaking within and beyond the SM are expected to dominate the scene of particle-physics 
phenomenology in the present and next decades. 

One of the most theoretically appealing realizations of the Higgs mechanism for mass 
generation is provided by Supersymmetry (SUSY). The minimal super symmetric exten- 
sion of the SM (MSSM) has a number of interesting field-theoretic and phenomenological 
properties, if SUSY is softly broken such that superparticles acquire masses not greatly ex- 
ceeding 1 TeV. Specifically, within the MSSM, the gauge hierarchy can be made technically 
natural Unlike the SM, the MSSM exhibits quantitatively reliable gauge-coupling uni- 
fication at the energy scale of the order of 10^^ GeV jl]. Furthermore, the MSSM provides 
a successful mechanism for cosmological baryogenesis via a strongly first-order electroweak 
phase transition jSHEl; and provides viable candidates for cold dark matter [71IH]. 

The MSSM makes a crucial and definite prediction for future high-energy experiments, 
that can be directly tested at the Tevatron and/or the LHC. It guarantees the existence 
of (at least) one light neutral Higgs boson with mass less than about 135 GeV 0. This 
rather strict upper bound on the lightest Higgs boson mass is in accord with global analyses 
of the electroweak precision data, which point towards a relatively light SM Higgs boson, 
with Mfjg^ ~ 211 GeV at the 95 % confidence level (TU]. Furthermore, because of the 
decoupling properties of heavy superpartners, the MSSM predictions for the electroweak 
precision observables can easily be made consistent with all the experimental data [TT] . 

Recently, a new important phenomenological feature of the MSSM Higgs sector has 
been observed. It has been realized that loop effects mediated dominantly by third- 
generation squarks may lead to sizeable violations of the tree-level CP invariance of the 
MSSM Higgs potential, giving rise to significant Higgs scalar-pseudoscalar transitions |12j . 
in particular. As a consequence, the three neutral Higgs mass eigenstates -^1,2,3, labelled 
in order of increasing mass such that Mh^ < < Mh^, have no definite CP parities, 
but become mixtures of CP-even and CP-odd states. Much work has been devoted to 
studying in greater detail this radiative Higgs-sector CP violation in the framework of the 
MSSM |Il[Tl[TnilISllIZllIHllIll2niinil2Sll221- In the MSSM with explicit CP violation, the 
couplings of the Higgs bosons to the SM gauge bosons and fermions, to their supersymmet- 
ric partners and to the Higgs bosons themselves may be considerably modified from those 
predicted in the CP-conserving case. Consequently, radiative CP violation in the MSSM 
Higgs sector can significantly affect the production rates and decay branching fractions of 
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the Higgs bosons. In particular, the drastic modification of the couphngs of the Z boson to 
the two fighter Higgs bosons B.\ and migfit enable a relatively ligfit Higgs boson witfi a 
mass Mii^ even less tfian about 70 GeV to fiave escaped detection at LEP 2 [THll22j . Tfie 
upgraded Tevatron collider and tfie LHC will be able to cover a large fraction of tfie MSSM 
parameter space, including tfie cfiallenging regions witfi a ligfit Higgs boson witfiout definite 
CP parity [SllESESI • Furtfiermore, complementary and more accurate explorations of tfie 
CP-noninvariant MSSM Higgs sector can be carried out using fiigfi-luminosity e^e~ [22] 
and/or 77 j2Zj colliders. In addition, a complete determination of tfie CP properties of tfie 
neutral Higgs bosons is possible at muon colliders by exploiting polarized muon beams |28j . 

It is obvious tfiat a systematic study of Higgs pfienomenology in tfie MSSM witfi 
explicit CP violation would be greatly facilitated by an appropriate computational tool. 
For tfiis purpose, we fiave developed tfie Fortran program CPsuperH, a new self-contained 
computational package, wfiicfi calculates tfie mass spectrum and tfie decay widtfis of tfie 
neutral and cfiarged Higgs bosons in tfie MSSM witfi explicit CP violation^^. It calcu- 
lates tfie neutral Higgs-boson masses Mh^ ^ .^ and tfie corresponding 3x3 Higgs-boson 
mixing matrix O, employing tfie renormalizat ion-group- (RG-)improved diagrammatic ap- 
proacfi of |2S|- We include tfie leading two-loop QCD logaritfimic corrections as well as tfie 
leading two-loop logaritfimic corrections induced by tfie top- and bottom-quark Yukawa 
couplings [201 • We also include tfie leading one- loop logaritfimic corrections due to gaugino 
and fiiggsino quantum effects [21]. Moreover, we implement tfie potentially large two-loop 
non-logaritfimic corrections originating from one-loop tfiresfiold effects on tfie top- and 
bottom-quark Yukawa couplings, associated witfi tfie decoupling of tfie tfiird-generation 
squarks [S21CS1E]- Finally, tfie RG-improved diagrammatic calculation takes account of 
mass sfiifts determined by tfie positions of tfie poles in tfie corresponding Higgs-boson prop- 
agators. Tfiese Higgs-boson pole-mass sfiifts for tfie ligfitest Higgs boson are small, of tfie 
order of a few GeV. However, tfie mass sfiifts for tfie fieavier Higgs bosons, and ifa, 
can be mucfi larger and of tfie order of several tens of GeV [22]; especially if tfieir masses 
Mii^ and M^.^ fiappen to be close to tfiresfiolds for tfie on-sfiell production of squark pairs. 
Finally, we note tfiat tfie computation of all tfie Higgs-boson decay widtfis by tfie code 
CPsuperH relies on tfie extensive analytic results for tfie decay widtfis presented in [TH] . 

After tfiis introductory discussion, tfie next section summarizes all tfie topics of Higgs 
pfienomenology tfiat can be studied witfi tfie code CPsuperH. In Section [21 we describe tfie 
execution procedure of tfie Fortran code and present examples of input and output files 
from a test run. Finally, we summarize tfie essential features of tfie code and provide an 
outlook for furtfier developments of CPsuperH in Section |3] 

'''We note in passing that a Fortran code called HDECAY has already been developed for the calculation 
of Higgs boson decays in the CP-invariant version of the MSSM [21| . 
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2 Higgs Phenomenology in the MSSM with Exphcit 
CP Violation 



In the presence of nontrivial CP-violating phases for the higgsino mass parameter /i and 
the soft SUSY-breaking parameters in the MSSM, the couphngs of the Higgs bosons to the 
gauge bosons, fermions and sfermions, as well as those to the Higgs bosons themselves, are 
strongly modified. In order to investigate these modifications and their phenomenological 
implications, we begin by stating our conventions for the mixing matrices of the neutral 
Higgs bosons and SUSY particles, and then we present all the relevant Higgs interactions 
with the MSSM particles to be used subsequently for calculating the masses, total decay 
widths and decay branching fractions of the neutral and charged Higgs bosons. 



2.1 Conventions 

In this subsection, we give our conventions for the mixing matrices of neutral Higgs bosons, 
charginos, neutralinos and third-generation sfermions. 

• Neutral Higgs bosons: In the presence of nontrivial CP-violating phases of the soft 
supersymmetry-breaking parameters, most relevantly in the third-generation sfermion 
sector, the three neutral Higgs bosons all mix together via radiative corrections: 

{(P^,4>2,a)l = 0^,{H,,H2,Hs)J, (1) 

where O^MjjO = diag (M^^, M^^, M^J with Mh, < Mh, < Mh,. We refer 
to [iniEHl for the details of the calculations of the mass-squared matrix M^h-, the 
diagonalization matrix O and the pole masses of the Higgs bosons. 

• Charginos: In SUSY theories, the spin-1/2 partners of the gauge bosons and 
the charged Higgs bosons, and if^, mix to form chargino mass eigenstates. We 
adopt the convention Hj^^^^ = H^^) ■, where the subscripts 1 and 2 are associated with 
the Higgs supermultiplets leading to the tree-level mass generation of the down- and 
up-type quarks, respectively. The chargino mass matrix in the {W~,H~) basis 

Ma V2Mw cp \ 

(2) 

V 2Miy SfS J 

is diagonalized by two different unitary matrices CrJ^cCl = diag{m^±, fn<^±}, where 
m-± < fn^±. The chargino mixing matrices {Ci)ia and {CR)ia relate the electroweak 
eigenstates to the mass eigenstates, via 

XaR = iCR)lXm, XaR = (3) 
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M 



N 



(4) 



We use the following abbreviations throughout this paper: sp = sin (3, cp = cos 
= tan/3, S2/3 = sin 2/5, C2/3 = cos 2/3, sw = sin 9w, cw = cos 9w, etc. 

Neutralinos : The neutralino mass matrix in the {B, W^, H^, H2) basis is given by 

f Ml -Mzcpsw Mzspsw \ 

M2 MzCpCw -MzSpCw 

-Mzcpsw Mzcpcw -/X 

V Mzspsw -Mzspcw -fi y 

This neutralino mass matrix is diagonalized by a unitary matrix A^: N*M.nN''' = 
diag (m^o, m^o, m^o, m^o) with m^o < m^o < m^o < m^o. The neutralino mixing 
matrix Nia relates the electroweak eigenstates to the mass eigenstates via 

(5,ly^^?,^°)I = iV*,(x?,x^,xg,xD^ (5) 

Stops, sbottoms, staus and tau sneutrino: The stop and sbottom mass matrices may 
conveniently be written in the {qi, qn) basis as 



w) 



hlv,{Al-f,R,)/V2 



w 



with q = t,b, R = U,D, Tl = -T^ = 1/2, Qt = 2/3, Qb 



-1/3, Vb 



(6) 



V2, 



Rb = tan/3 = V2/V1, Rt = cot /3, and hq is the Yukawa coupling of the quark q. On 
the other hand, the stau mass matrix is written in the {fi, fji) basis as 



Mt 



Ml +ml + C2pMl {si, - 1/2) K<A*^ - /itan/3)/V2 



hrVi{Ar - /i*tan/3)/y2 



^4, 



C2/3MI siy 



(7) 



and the mass of the tau sneutrino is simply m^^ = yM? + |c2/3M|, as it has 
no right-handed counterpart in the MSSM. The 2x2 sfermion mass matrix for 
f = t,b and r is diagonalized by a unitary matrix U-^: U^'^ Ai'i = diag( 



h 12 



with < m|^. The mixing matrix relates the electroweak eigenstates Jl^r to 
the mass eigenstates /i,2, via 

{fLj~R)l = uLifuM. (8) 

We parameterize the mixing matrices as follows: 

cos 6 
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sin^^-e+*'^/ 



-smOje-^'^'f 
cos 9? 



(9) 



and we calculate numerically the mixing angle and phase 0^ in the ranges between 
— 7r/2 and 7r/2, so that cos6j > and cos0j > 0. 
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2.2 Higgs— Boson Interactions 



In this subsection, we list all the Higgs interactions with gauge bosons, SM fermions, 
squarks, sleptons, charginos, and neutralinos. We also present all the trilinear and quartic 
Higgs-boson self-couplings. 

• Interactions of Higgs bosons with gauge bosons: The interactions of the Higgs bosons 
with the gauge bosons Z and are described by the three interaction Lagrangians: 

Luvv = gMw (w;W-^ + ^^M^^) T.9n^yyH,, (10) 

ChHZ = -^Y.9n^n,zZ'\H,l'd,H,), (11) 

Chh±w^ = -| 5: 9,^„^^. W->^{H,id, + h.c. , (12) 

i 

where g = e/ sin6'vy is the SU(2)l gauge-coupling constant, and the couplings g^-wy 
(7„ „ ^ and q ^ are given in terms of the neutral Higgs-boson mixing matrix O 
by (note that det(O) = ±1 for any orthogonal matrix O): 

9h,h^z = sign[det(0)] e^j-fc^^^^^ , 
g„^„+^- = Cf3 O^^, - 0^,i - iOai , (13) 

leading to the following sum rules: 

3 

Efif^ =1 and + \q . P = 1 for each 2 . (14) 

1=1 

• Higgs-quark-antiquark and Higgs-lepton-antilepton interactions: The effective La- 
grangian governing the interactions of the neutral Higgs bosons with quarks and 
charged leptons is 



^Hjf = - E E / (4/7 + </775) / • (15) 

f=u,d,l ^ *=1 



At the tree level, {g^, g^) = {O^^i/ c^, -Oai tan (3) and {g^ , g^) = {O^^i/ sp, -Oai cot (3) 
for / = (/, d) and f = u, respectively. In the case of third-generation quarks, the pro- 
gram CPsuperH computes the finite threshold corrections induced by the exchanges 
of gluinos and charginos. As described in Appendix A, we include the all-orders 
resummation j221Elll22] of the leading powers of tan/3, as required for a meaning- 
ful perturbative expansion. Correspondingly, in the presence of CP violation, the 
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effective couplings of the charged Higgs boson to quarks and leptons in the weak- 
interaction basis are described by the interaction Lagrangian: 

V^^^{/T,/i)={M),(-,0 

(16) 

where Pl/r = (1 =F 75)/2- At the tree level, = cot P and = tan/3 with 
mj, = 0. The loop-induced threshold corrections to the g^+y^^ couplings are presented 
in Appendix A. 

• Higgs-sfermion-sfermion interactions: The Higgs-sfermion-sfermion interactions can 
be written in terms of the sfermion mass eigenstates as 

^H// = - E 9uJ^fSH^fU^^ (17) 

f=u,d 

where 




with a = (01, 02, a) = (1,2,3), p,-f = L,R, i = {Hi,H2,H^) = (1,2,3) and 
j,k = 1,2. Likewise, the charged Higgs-boson interactions with up- and down-type 
sfermions are given by 

^H±ff' = ^9H-^f*rSH^ /; ~f'k) + h.c, (18) 

where 




The expressions for the couplings T^^''^ and T^^^*^' to third-generation sfermions are 
presented in Appendix B. As shown in [21], our iterative treatment of the threshold 
corrections that are enhanced at large tan (3 ensures that the corresponding corrections 
to the Higgs-boson couplings to squarks are also resummed correspondingly. 

• Interactions of neutral Higgs bosons and charginos: These are described by the fol- 
lowing Lagrangian: 

V i,j,k 

4,x+x7 = ^{[(C'K)a(C'i)*2G'r + (C'«).2(Ci)*iGf^] + [z-jr} , 
9^^^ = I {[{CpMCL)*2Gt + {CRh{CL)*iGt] - [^ - J?} , (19) 
where Cf.' = {O^^k - ispOak), Gt = (O^^k - icpOak), i,j = 1, 2, and A; = 1 - 3. 
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Interactions of neutral Higgs bosons and neutralinos: These are described by the fol- 
lowing Lagrangian: 



j)] 



--^m[{N;, - twN*,){N*,Gt - N*,Gt) + (z ^ j)] , (20) 



where i,j = 1-4 for the four neutralino states and k = 1-3 for the three neutral Higgs 
bosons. 

Interactions of charged Higgs bosons, charginos and neutralinos: These are described 
by the following Lagrangian: 



3h+x°x- 



J2 



C/3 



V2N,^iGR)*, + (iV,2 + twNa)iGR)*^] } , 



J2 



C/3 



V2N,,iGn)*, + {N,,2 + twNn)iGj,)*,] } . 



(21) 



Trilinear and quartic Higgs-boson self-couplings P^l^ : The effective trilinear and 

(22) 



quartic Higgs self-couplings can be cast into the form 

3 3 



3H 



i>j>k=l i=l 



^4H = E 9h,hjHi^h, HiHjHkHi + E 9h.h H+H- HiHjH+H' 



where 



i>j>k>l=l 

9, 



i>j=l 



-H+H- ^ 



E {OaiOpjO-yk} galS'y , 9 



HiH+H- 



E 



H+H- ' 



(23) 



(24) 



Q</3<7=1 
3 

E {OaiOi3j0^kOsi} QaP^S , 
a</3<7<<5=l 

3 

9af3H+H- ■ 

a<f3=l 



a=l 



(25) 
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In the above equations ()24|1 and (j2SI), the expressions within the curly brackets {■ ■ ■} 
need to be symmetrized with respect to the indices i,j,k,l and divided by the cor- 
responding symmetry factors in cases where two or more indices are the same. For 
example, {OaiOisjO^k} can explicitly be evaluated as follows: 



{OaiOfjjO^k} = 




+OakO^,0^j + O^kO^jO^i) , (26) 

with Ns = 6 when i = j = k, Ns = 1 when k) = (3, 2, 1), and Ns = 2 in all 
the other cases. We present the couplings gaP'y, QaH+H-^ QaP'^s^ and QajdH+H- of the 
Higgs weak eigenstates in Appendix C. 



2.3 Neutral and Charged Higgs Boson Decays 

In this subsection, we calculate all the two-body decay widths of the Higgs bosons. We 
consider the decays of the Higgs bosons into pairs of leptons, quarks, charginos, neutralinos, 
massive gauge bosons, Higgs bosons, squarks, sleptons, photons and gluons as well as into 
a massive gauge boson and a Higgs boson. For the decay modes involving more than one 
massive gauge boson, three-body decays are also considered [TH] . 



f f: First, let us consider the decays into a pair of fermions. Without loss 



of generality, the Lagrangian describing the interactions of the Higgs bosons with two 
fermions can be written as 



C 



T.9frH^fk{ 



,75)/; +h.c. 



(27) 



where /*^'^ stands for a lepton, a quark, a chargino, or a neutralino, and the tree-level 
couplings Qf, Qff and g^'^ are given in Tabled In terms of these generic couplings, 
the width for a decay into two Dirac fermions is given by 



D 



N, 



in 



■i^k){\9 



S\2 



\9 



2y/KjKk{\g 



S\2 



\9 



P|2^ 



(2^ 



where k,- 



c 



mjjMjj and X{l,x,y) = (1 - x 
-iPk\9s? + Ifi-pP) when Kj = Kk 



uY — 4:xy. We note that Yd becomes 

: K, where (3^ = Vl ~ 4k. The colour 
factor Nc = 3 for quarks and 1 for leptons, charginos, and neutralinos. The decay 
widths into two Majorana fermions are given by 

/ 4 \ „ 



M 



1 + 6. 



D 



(29) 
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where 6jk = 1 for identical Majorana fermions. 



Table 1: The couplings gj, gjf and g^'^ in Eq. \21\^ at the tree level. 



Decay Mode 


S' P 

9f 9 9 


Hi -> l+l- 
Hi dd 
Hi uu 

Hi XjXk 


S 0^,^/Sp {Cfs/Sf,)Oa^ 

4x|x.- 4xtx.- 


Decay Mode 


9ff' 9^ 9^ 


H+ l+u 
H+ ud 
H^ -> xU'i 


-Jf^ [lltp + {mdlmu)tp\l2 z[l/t^-(md/m„)ts]/2 
9l^ 9h+x°x- ^h+xHx- 



For the couplings S'^^f,, Qh^ii and g^+tf,, the finite loop-induced threshold correc- 
tions due to the exchanges of gluinos and charginos can be included by taking 
IFLAG_H(10)=0 (the default setting) in the code CPsuperH, as explained in Sec. El 

For Hi qq, the leading-order QCD correction is taken into account by applying the 

enhancement factor Kf = 1 -|- 5.67 ^"'^ to the decay width given above. We take 
the running fermion masses at the scale mf"^*^ as reference values. The effect on the 
couplings of the running of the quark masses from the top-quark pole mass scale to 
the Higgs-boson mass scale is also considered in calculating the corresponding decay 
widths, as 

qq) = Kl f ^f^f Vr^(i/. qq) , (30) 

where mq(mf°''') is used in gf, but mq{mHi) is used in kj, when calculating Tq. 
Likewise, running q and q' quark masses are used when computing T{H^ gg'), 
while the dominant one-loop QCD corrections have been included by factors very 
similar to Kf. Finite quark- mass and higher-order QCD effects will depend, to some 
extent, on the CP-violating parameters of the MSSM, and require an independent 
study. In the present version of the code, we only include the leading-order QCD 
effects which remain unaffected by CP violation. Also, we do not include flavour- 
violating decays of the neutral Higgs bosons. We plan to implement such refinements 
in future versions of CPsuperH. 
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Hj VV: The width for decay into two massive gauge bosons is given by 



Gpg^ Ml 5v 
167271 



Ay(l-4«:,y + 12«:?y), (31) 



where Gf/\/^ = g'^/8M^, Kiv = M^/Mfj^, /3iv = y/l - ^i^iv, and = 2 and 5z = 
M^/{cwMzY = 1. The three-body decay width r{Hi VV*) is also calculated, 
using 



r(H VV*) = ^P3%vv^hMyv' f(V^-i)^ ^^ev\"\uj,,x, l)[A(.7„x, 1) + 12x] 



16v^7r h cuf7r[(x-l)2 + e2,] 

(32) 

where uji = MfjjMy = l/niv, Syv* = 2, and ey = Ty/My. We note that 
\^/\u:i, 1, l)[A(u;„ 1, 1) + 12]/a;f = Ay(l - 4«:.y + 12^%) . 

Hi HjZ and HjW^: The decay width of a heavier Higgs boson into a 

lighter Higgs boson and a massive gauge boson is given by 



Sv^tt'"'-' ,.,,.„kv), (33) 

where {MH,g) = {MH,,g,^„^,) or {Mh+, g^^^^^_), = MljM^, and Ky = 
M^/Mjj. The three-body decay widths r(if/^ HjZ*) and T{H+ HjW+*) are 
given [IH] by: 

r(ff. H,Z') . ''J^A^ d. ff <^:r'-1, , (34) 

^ ^ ^ 8V27r k wf7r[(x - 1)2 + e|] ' ^ ' 

with cijj = M\JM\^ and similarly 



8V27r cui7r[(a; - 1)2 + e^] ' 

with ujj = MjjjM'^ and = Mfj^/W^. 

~ — ~ — * 

/Jj — s> HjHk, Hi fjfki and H^ fjf'k'- The decay widths into two scalar particles 

can be written as 

r = Nf ^A^/2(l,K„Kfe), (36) 

where* {Nf,Q) = (l+^jfc, ^h^h^hJ, i^c^gnj*/,) or {Nc, gn+f- f'J, and k^- = Ml^ jjM^.. 



^3 ^J3 

^We note that the couplmgs gH.HjH^^ ^^'^ defined via the Lagrangian The HiHjHj vertices for 



i > j therefore contain an extra factor of 2. 



11 



Hi 77: The amplitude for the decay process Hi — > 77 can be written as 

M,,H. = -^{sUMhJ (etx ■ eL) - P^iMnJj^ielelhh)^ , (37) 

where ki^2 are the momenta of the two photons and ei^2 the wave vectors of the 
corresponding photons, eij_ = — 2A;f(A;2 ■ ^i)/^Hi^ ^2± = ^2 2/c2(^i • ^2)/M]j_ 
and (eie2/ci/c2) = e/^i/po- £2^1/^2 • The scalar and pseudoscalar form factors, retaining 
only the dominant loop contributions from the third-generation (s)fermions, and 
charged Higgs bosons, are given by 



P^{Mh^) = 2 Y: NcQ}gfg^jf^F,j{T,f), (38) 

where Ti^ = M^./4m^, Nq = 3 for (s) quarks and Nc = 1 for staus and charginos, 
respectively. 

The form factors Fgf, Fpf, Fq, and Fi can be expressed in terms of a so-called scaling 
function /(r), by 

Fsfir) = r-Ml + (l-r-i)/(r)], F,j{r) = f{r) , (39) 
Fo(r) = r-M-l + r-V(r)], F^ir) = 2 + 3r-' + 3r-\2 - r-')f{r) , 

where /(r) stands for the integrated function 

1 /•! dw f arcsin^(A/r) : t<1, 

f(r) = — Jiin[l-ATy(l-y)] = { ,r . /-^ n2 ~ (40) 

2^0 y ^ I -i [ln( ^+:^ ) -ztt] : r>l. ^ 

It is clear that imaginary parts of the form factors appear for Higgs-boson masses 
greater than twice the mass of the charged particle running in the loop, i.e., r > 1. 
In the limit r ^ 0, F,/(0) = 2/3, Fpf{0) = 1, Fo(0) = 1/3, and Fi(0) = 7. Finally, 
the decay width is given by 

/Vf3 ^2 

FiH, 77) = TT?:^ + iPHMnjf] . (41) 



2567r'^f^ - 

The QCD correction to the width T{Hi 77) is included in the large loop-mass limit 
by multiplying the rescaling and factors to the quark and squark contributions 
to SJ, respectively. The rescaling factors in the large loop-mass limit are given by [SS] 



J2 = \- ^ , Jl = l+ ' • (42) 

TT 37r 
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Hi gg: The amplitude for the decay process Hi ^ gg {i = 1, 2, 3) can be written 



as 



M,,H. = - {^f(M^J (e^x ■ - P!{Mn:)j^{e\elhk,)] , (43) 

where a and 6 (a, 6 = 1 to 8) are indices of the eight SU(3) generators in the adjoint 
representation, and ki^2 and ei,2 are the four-momenta and wave vectors of the two 
gluons, respectively. The scalar and pseudoscalar form factors, retaining only the 
dominant contributions from third-generation (s) quarks, are given by 



/=^'* ^ fj=ii,i2M~b2 fj 

P!{Mnd = T.9f9^ff-^F,j{nf). (44) 
f=b,t '"-f 

The decay width of the process Hi gg is then given by 

]U3 2 

- 99) = tt^tH ISfiMnjf + Al iPfiMnjf] , (45) 

where Kfj ^ are QCD loop enhancement factors that include the leading-order QCD 
corrections. In the heavy-quark limit, the factors Kfj ^ are given by 

TT V 4 



TT V 4 6 



- Nf), (46) 



where Np is the number of quark flavours lighter than the Hi boson. Away from 
isolated regions of the parameter space [33] , the above Kf^ j^^ factors capture the main 
bulk of the NLO corrections with an accuracy at the 10% level. 



3 The Structure of CPsuperH 

The program CPsuperH is self-contained with all necessary subroutines included. The 
Fortran code CPsuperH uses three input arrays for reading the input parameters and five 
output arrays for generating the Higgs couplings, decay widths and branching fractions. 
The three input arrays are named 

SMPARA_H(IP), SSPARA_H(IP), IFLAG_H(NFLAG) . 

^We ignore the small difference between the i^— factors of the quark and squark loop contributions to 
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Among the five output arrays, the arrays for generating the Higgs couphngs are named 



NHC_H(NC,IH), SHC_H(NC), CHC_H(NC), 

and the other two output arrays for the decay widths and branching fractions are named 

GAMBRN ( IM , IWB , IH) , GAMBRC ( IM , IWB) . 

The code CPsuperH contains also arrays for the masses and mixing matrices of the Higgs 
bosons, stops, sbottoms, staus, charginos and neutrahnos as explained below. 

3.1 Input Arrays 

In this subsection, we describe the details of the input arrays. 

• SMPARA_H(IP) : This is the array for the SM input parameters. In the current version, 
we are dealing with 15 inputs as shown in Table |2l but these can easily be extended 
by changing NSMIN in cpsuperh.f . This array is filled from the file run . 



Table 2: The contents o/ SMPARA_H(IP) . 



IP 


Parameter 


IP 


Parameter 


IP 


Parameter 


IP 


Parameter 


1 

2 
3 
4 

5 


asiMz) 

Mz 
sin^ 6w 

me 


6 
7 
8 
9 
10 


nij. 

1 polc\ 

md[mt ) 


11 
12 
13 
14 
15 


pole 

Tz 


16 
17 
18 
19 
20 





• SSPARA_H(IP) : This array is for the SUSY input parameters. In the current version, 
we are dealing with 21 inputs as shown in Table El but these can easily be extended 
as well by changing NSSIN in cpsuperh.f. This array is also filled from the file run . 

• IFLAG_H(NFLAG) : This NFLAG-dimensional array controls CPsuperH. This fiag array 
is used for printing options, calculating options, integer input parameters, error mes- 
sages, etc. The default value for every fiag is zero. This array also can be filled from 
the file run . Only a part of IFLAG_H is being used presently by the code: 

— IFLAG_H(1)=1: Print out the input parameters. 

— IFLAG_H(2)=1: Print out the masses and mixing matrix of the Higgs bosons. 
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Table 3: The contents o/ SSPARA_H(IP) . 



IP 


Parameter 


IP 


Parameter 


IP 


Parameter 


IP 


Parameter 


IP 


Parameter 


1 

2 
3 
4 
5 


tan / 
MP° 

% 
\Mi 


3 

Ic 

b 


6 
7 
8 
9 
10 


$1 
IM2I 
$2 

IM3I 

$3 


11 

12 
13 
14 
15 




16 
17 
18 
19 
20 




At 
A, 

Ar 




21 
22 
23 
24 
25 





— IFLAG_H(3)=1: Print out the masses and mixing matrices of the stops, sbottoms, 
tau sneutrino and staus. 

— IFLAG_H(4)=1: Print out the masses and mixing matrices of the charginos and 
neutrahnos. 

— IFLAG_H(5)=IX: Print out the Higgs-boson couphngs. The couphngs oi Hi, H2, 
H3 and to two particles will be printed for IX =1, 2, 3, and 4, respectively, 
and the Higgs-boson self-couplings will be printed for IX=5. All these couplings 
can be printed out altogether by taking IX=6. 

— IFLAG_H(6)=IX: Print out the decay widths and branching ratios. The decay 
widths and branching ratios of Hi, H2, H3, and H^ will be printed for IX =1, 
2, 3, and 4, respectively. IX=5 is for printing out all the decay widths and 
branching ratios of the neutral and charged Higgs bosons. 

— IFLAG_H(10)=1: Do not include the finite threshold corrections to the top- and 
bottom-quark Yukawa couplings due to the exchanges of gluinos and charginos. 

— IFLAG_H(11)=1: Use the effective potential masses for Higgs bosons instead of 
their pole masses. 

— IFLAG_H(20)= ISMN. The index ISMN is used for GAMBRN, which is an array for the 
neutral Higgs decay widths, and its default value is ISMN = 50. In general, (ISMN- 
1) is the maximal number of different decay modes of the neutral Higgs bosons 
into SM particles. The value of ISMN may be changed in order to incorporate 
additional rare decay modes. The index ISMN is reserved for the subtotal decay 
width and branching fraction of the decays into the SM particles in the output 
array GAMBRN (see below). 

— IFLAG_H(21)=ISUSYN. Similarly to ISMN, the index ISUSYN is used for GAMBRN, 
and its default value is ISUSYN = 50, with (ISUSYN-1) being the maximal num- 
ber of different decay modes of the neutral Higgs bosons into SUSY particles. 
The index ISMN+ISUSYN is reserved for the subtotal decay width and branching 
fraction of the decays into the SUSY particles, while the index ISMN+ISUSYN+1 
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is used for the total decay width, considering decays into both the SM and SUSY 
particles in the output array GAMBRN (see below). 

- IFLAG_H(22)= ISMC. The index ISMC is used for GAMBRC, which is an array for 
the charged-Higgs decay width, and its default value is ISMC = 25. In general, 
(lSMC-1) is the maximal number of different decay modes of the charged Higgs 
bosons into SM particles. The index ISMC is reserved for the subtotal decay 
width and branching fraction of the decays into the SM particles in the output 
array GAMBRC (see below). 

- IFLAG_H(23)=ISUSYC. Similarly to ISMC, the index ISUSYC for GAMBRC, with 
(lSUSYC-1) being the maximal number of different decay modes of the charged 
Higgs bosons into SUSY particles. The index ISMC+ISUSYC is reserved for the 
subtotal decay width and branching fraction of the decays into the SUSY parti- 
cles, while the index ISMC+ISUSYC+1 is used for the total decay width considering 
decays both into the SM and SUSY particles in the output array GAMBRC (see 
below) . 

In Appendix D, we list all the parameter common blocks filled or calculated from 
SMPARA_H and SSPARA_H. 

3.2 Output Arrays 

In this subsection, we give detailed descriptions of the output arrays. Some of the entries 
of IFLAG_H are reserved for various error messages. This feature might be helpful when 
using CPsuperH to scan many parameter points: 

• IFLAG_H(50)=1: This is an error message that appears when a stop or sbottom 
squared mass is negative. 

• IFLAG_H(51)=1: This is an error message that appears when the Higgs-boson mass 
matrix contains a complex or negative eigenvalue. 

• IFLAG_H(52)=1: This is an error message that appears when the diagonalization of 
the Higgs mass matrix is not successful. 

• IFLAG_H(53)=1: This is a warning message that appears when the second-step im- 
provement in the calculations of the pole masses is needed. 

• IFLAG_H(54)=1: This is an error message that appears when the iteration resumming 
the threshold corrections is not convergent. 

• IFLAG_H(55)=1: This is an error message that appears when the Yukawa coupling 
has a non-perturbative value: \ht\ or \hi,\ > 2. 
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IFLAG_H(56)=1: This is an error message that appears when a tau sneutrino or a 
stau squared mass is negative. 

The main numerical output is stored in the following arrays: 

NHC_H(NC,IH): This is an array for the IH-th neutral Higgs boson {Hja) couplings to 
two particles with index NC. Currently, this array is filled up to NC = 93 as shown in 
Table i| 

SHC_H(NC): This array is for the self-couplings of Higgs bosons. Currently, this array 
is filled up to NC = 35, as shown in Table |31 

CHC_H(NC) : This array is for the couplings of the charged Higgs boson to two particles. 
Currently, this array is filled up to NC = 48, as shown in Table IHl 

GAMBRN(IM,IWB,IH): This output array is for the decay width in GeV (IWB=1) and 
branching fraction (IWB=2,3) of the decay mode specified by the index IM of the 
neutral Higgs bosons Hj^. The value IWB=2 is for the branching fraction taking 
into account the decays only into SM particles, and IWB=3 for that taking account 
both the SM and SUSY decays. By default, the code takes ISMN = ISUSYN = 50. 
All the decay modes considered are listed in Table U\ for specific IM. In particular, 
GAMBRN(IM=ISMN+ISUSYN+1 , IWB=1 , IH) is the total decay width of the neutral Higgs 
boson Hm and GAMBRN(IM=ISMN, IWB=1 , IH) and GAMBRN(IM=ISMN+ISUSYN, IWB=1 , IH) 
are the subtotal decay widths into SM particles and into SUSY particles, respectively. 
Therefore, we have the following relations for the branching fractions IWB=2,3 : 

GAMBRN(ISMN, 2, IH) = GAMBRN(ISMN + ISUSYN + 1, 3, IH) = 1 

and 

GAMBRN(ISMN, 3, IH) < 1 , GAMBRN(ISMN + ISUSYN + 1, 2, IH) > 1. 

GAMBRCdM, IWB) : This array is for the decay width in GeV (IWB=1) and branching 
fraction (IWB=2,3) of the decay mode number IM of the charged Higgs boson. The 
convention for IWB is the same as that for GAMBRN. In the code, ISMC = ISUSYC = 
25 is taken. The decay modes considered are shown in Table |H1 In particular, 
GAMBRC(IM=ISMC+ISUSYC+1 , IWB=1) is the total decay width of the charged Higgs bo- 
son and GAMBRC(IM=ISMC, IWB=1) and GAMBRC(IM=ISMC+ISUSYC, IWB=1) are subtotal 
decay widths into SM and into SUSY particles, respectively. Similarly to the case of 
the neutral Higgs bosons, we have the relations 

GAMBRC(ISMC, 2) = GAMBRC(ISMC + ISUSYC + 1, 3) = 1 
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and 

GAMBRC(ISMC, 3) < 1 , GAMBRC(ISMC + ISUSYC + 1, 2) > 1. 

• The code CPsuperH contains output arrays for the masses and mixing matrices of the 
neutral Higgs bosons, the sfermions, the charginos, and the neutrahnos, named as 
follows: 

— HMASS_H(3): The masses of the three neutral Higgs bosons, Mh^- 

— 0MIX_H(3,3): The 3x3 Higgs mixing matrix, O^i- 
~ STMASS_H(2): The masses of the stops, m^^. 

— STMIX_H(2,2): The mixing matrix of the stops, U^^. 

— SBMASS_H(2): The masses of the sbottoms, m^.. 

— SBMIX_H(2,2): The mixing matrix of the sbottoms, t/^j. 

— STAUMASS_H(2): The masses of the staus, nif^. 

— STAUMIX_H(2,2): The mixing matrix of the staus, U^^. 

— SNU3MASS_H: The mass of the tau sneutrino, mc,^. 

— MC_H(2): The masses of the charginos, m-±. 

— UL_H(2,2): The mixing matrix of the left-handed charginos, (Cl)^^. 

— UR_H(2,2): The mixing matrix of the right-handed charginos, {CRjia- 

— MN_H(4): The masses of the neutralinos, m^o. 

— N_H(4,4): The mixing matrix of the neutralinos, Ni^- 

3.3 How to Run CPsuperH 

The package CPsuperH consists of two text, five Fortran, and three shell-script files. The 
main features of the files follows: 

• Text files : 

— The file ARRAY shows all the arrays described in the previous two subsections. 

— The file COMMON lists the parameter common blocks, as described in Appendix 
D. 

• Fortran files : 

— cpsuperh . f fills all the arrays in ARRAY from the shell-script file run by calling 
the following four Fortran files. 
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— f illpara.f fills the common blocks in COMMDN from SMPARA_H. and SSPARA_H. 

— f illhiggs . f fills the arrays for the masses and the mixing matrix of the neutral 
Higgs bosons, HMASS_H and OMIX_H. 

— f illcoupl . f fills the arrays for the masses and the mixing matrices of the stops, 
the sbottoms, the charginos, and the neutralinos as well as the couplings arrays 
NHC_H, SHC_H, and CHC_H. 

— f illgambr.f fills the arrays GAMBRN and GAMBRC. 
• Shell-script files: 

— make lib creates the library file libcpsuperh. a from the four Fortran files of 
fillpara.f, fillhiggs.f, fillcoupl.f, and f illgambr . f . 

— compit creates the execution file cpsuperh . exe by compiling cpsuperh . f , linked 
with the library libcpsuperh. a. 

— run supplies cpsuperh . f with the input values for SMPARA_H and SSPARA_H and 
part of 1FLAG_H, and then shows the results by running cpsuperh.exe. The 
example presented in the present work is based on the so-called CPX scenario 
with = = ^Ar = ^3 = 90°, MsusY = 500 GeV, tan/3 = 5 and 
Mfjf = 300 GeV. Details may be found by inspecting the file. We note that, in 
the example, only IFLAG_H(l) = 1 is turned on initially. The user will have to 
edit run to choose new sets of parameters. The original version of run provides 
ample explanations of the various input parameters. 

It is straightforward to run the code CPsuperH. Type './makelib' and './compit' 
followed by './run': 

Run CPsuperH: ./makelib ./compit ./run 

and then one can see some outputs depending on the values of 1FLAG_H(1 — 6). 

In Appendix E, we show some sample outputs from a CPsuperH test run based on the 
CPX scenario. All the values for the input parameters used in the test run (IFLAG_H(1) = 1), 
the masses and mixing matrix of the neutral Higgs bosons (1FLAG_H(2) = 1), the masses 
and mixing matrices of the charginos and neutralinos (1FLAG_H(4) = 1), the lightest Higgs 
boson couplings (IFLAG_H(5) = 1), and the decay width and branching fractions of the 
lightest Higgs boson (IFLAG_H(6) = 1) are generated by taking the IFLAG_H values given in 
the parentheses. 

In order to check whether the code CPsuperH generates numerical outputs consistent 
with those provided by the code HDECAY j^lj in the CP-invariant case, we run the code 
CPsuperH in the 'maximal mixing' scenario: = v^^susy with the common SUSY 
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scale MsusY = Mq^ = Mfj^ = M^,^ = M^^ = M^^ = 1 TeV and = 100 GeV, by setting 
all the CP phases to zeros, but not including the threshold corrections. Fig. ^ shows the 
branching ratios and total decay widths of the MSSM Higgs bosons as functions of the Higgs 
boson masses. Although most of the parameter space presented in Fig. ^ is already ruled 
out by Higgs searches at LEP, we use it to present an effective comparison of the results of 
the CPsuperH code with those of the HDECAY code in the small-tan /3 regime. We find that 
the CPsuperH results are indeed consistent with those obtained by the code HDECAY. There 
are very few visible discrepancies, for example in B{Hi gg) for Mh^ ^ 80 GeV, which 
may be due in part to the improved calculation of the Higgs-boson mass spectrum and the 
mixing matrix in CPsuperH. 

We also show in Figs. |21and|nithe branching ratios, GAMBRN (IM , 2 , IH) , and total decay 
widths, GAMBRN (ISMN,1,IH) and GAMBRC(ISMC, 1 , IH) , found in the CPX scenario, which 
are consistent with the results previously reported in ^Hj- Finally, in Fig. |3]we illustrate 
the strong phase dependence of the Higgs boson decay widths into charginos and neutrali- 
nos. The left frame is for the lightest neutral Higgs boson, which for the given choice of 
parameters can only decay into two lightest neutralinos, i.e. GAMBRN(ISMN+ISUSYN,3, 1) = 
GAMBRN(ISMN+1 ,3, 1) . The right frame shows the branching ratios into superparticles of the 
heavier neutral Higgs bosons, GAMBRN(ISMN+ISUSYN,3,2) and GAMBRN(ISMN+ISUSYN,3,3), 
as well as the charged Higgs boson, GAMBRC(ISMC+ISUSYC,3), which only receive contribu- 
tions from chargino/neutralino final states. The results of this figure are in close agreement 
with ref. [H]. 

4 Summary and Outlook 

We have presented a detailed description of the Fortran code CPsuperH, a new computa- 
tional package for studying Higgs phenomenology in the MSSM with explicit CP violation. 
Based on recent RG-improved diagrammatic calculations j2Hj, the program CPsuperH com- 
putes the neutral and charged Higgs-boson masses as well as the 3x3 neutral Higgs-boson 
mixing matrix O in the presence of CP violation in the MSSM Higgs sector. Although 
the dominant one- and two-loop contributions to the Higgs-boson self-energies are incorpo- 
rated, there are still finite but subdominant two-loop contributions that may cause shifts 
of 3-4 GeV p^lHH] in the lightest Higgs-boson mass. These subdominant contributions 
can be estimated jHZI to be of comparable size with the dominant three-loop effects. Be- 
cause of the current lack of a detailed three-loop calculation, the subdominant two-loop 
contributions are not included in the present version of the code CPsuperH. 

In addition to the Higgs mass spectrum, the program CPsuperH computes all the 
couplings and the decay widths of the neutral Higgs bosons -^1,2,3 and the charged Higgs 
boson H~^, incorporating the most important quantum corrections ^lEl]. In particular. 
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the leading-order QCD corrections are included for the Higgs decays into photons, gluons 
and most hadronic channels, so that (with the possible exception of squark final states) the 
theoretical uncertainties in these decay modes are kept small. 

The Fortran code CPsuperH provides several options that can be selected in the in- 
put files for generating a variety of outputs for Higgs boson masses, Higgs decays and 
their respective branching fractions. Its structure offers the possibility of extending the 
list of input data to include flavour textures of off-diagonal trilinear and/or soft-squark 
masses, as these arise in some predictive schemes of soft SUSY breaking, e.g., in minimal 
supergravity models or in models with gauge and anomaly mediation. The code CPsuperH 
allows for straightforward extensions such as the additions of possible lepton- or quark- 
fiavour-violating decays of the Higgs bosons. Another possible extension is the inclusion 
of loop-induced absorptive phases, which allows to generate CP-odd rate asymmetries in 
Higgs boson decays (SHI- 

Radiative Higgs-sector CP violation in the MSSM has a wealth of implications for 
many different areas in particle-physics phenomenology. CP- violating phenomena mediated 
by Higgs-boson exchanges may manifest themselves in a number of low-energy observables 
such as the electron, neutron and muon electric dipole moments |H9 | l ^ l ^ l42j. They may 
also affect flavour-changing neutral-current processes and CP asymmetries involving K and 
B mesons jlSlESlIll]- Moreover, CP-violating Higgs effects may influence the annihilation 
rates of cosmic relics and hence the abundance of dark matter in the Universe |H] . Finally, 
an accurate determination of the Higgs spectrum in the presence of CP violation is crucial 
for testing the viability of electroweak baryogenesis in the MSSM. 

To conclude, the Fortran code CPsuperH can be used as a powerful and efficient com- 
putational tool in quantitatively understanding these various phenomenological subjects, 
which are inter-related within the framework of the MSSM with explicit CP violation. Even 
in the CP-conserving case, CPsuperH is unique in computing the neutral and charged Higgs- 
boson couplings and masses with equally high level of precision, and should be therefore a 
useful tool for the study of MSSM Higgs phenomenology at present and future colliders. 
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A Threshold Corrections 



The exchanges of gluinos and charginos give finite loop-induced threshold corrections to 
the Yukawa couplings hu^a, with the structure 

V2md 



hd 
hu 



1 

vcos(3 1 + {5hd/ hd) + ( A/id/ hd) tan /3 
\/2mu 1 



(A.i: 



wsin/9 1 + (Shu/hu) + (A/i„//i„) cot/3 ' 

modifying the couplings of the neutral Higgs-boson mass eigenstate to the scalar and 
pseudoscalar fermion bilinears as follows: 



9Hjd = 

+ Im 



4>ii 



+ Kd tan P ) cos /3 
Kd (tan^ + 1) 



Re 



O 



p 

dUidd 



-Re 
— Im 



1 + Kd tau/S 

tan/5 — 
1 + Kdiajif] 
i^d 

1 + tan j3 J cos /3 



Oai + Im 



1 + tan j3 J cos /3 



Kd tan /3 \ 0<^^ 



1 + Kd tan /3 / cos (3 



(A.2) 



Re 
+ Im 



1 + cot (3 J sin /3 
K„ (cot^/? + 1) 



Re 



1 + Ku cot /3 / sin /3 
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p 



-Re 
— Im 



1 + K„ cot [3 
cot j3 — Ku 



Oai + Im 



6l2 



1 + Ku cot /3 

1 + K„ cot P J sin /3 
In the above equations, we have used the abbreviation 

(Ahg/hg) 



Ku cot P \ O, 



1 + Ku cot P J sin /5 



P2« 



(A.3) 



(A.4) 



1 + i5hjhg) ' 

for q = u,d. Detailed expressions for the threshold contributions {Shg/hg) and [Ahg/hg) 
can be found in j22]- 

For the couplings of the charged Higgs boson to quarks, we have 

COt/5(l + Pu) - Ku 



9H+ud 



yn+ud 



1 + Ku cot P 

tan/3(l + p}) - K*, 
1 + K*^ tan P 
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(A.5) 
(A.6) 



The quantities Ku^d are given in ()A.4j) . whilst the quantities Ku^d and pu,d in (jA.Sj) and ()A.6j) 
are defined as follows: 



1 + {^hjh,) ' 1 + {6hjh,) 

The expressions for the additional threshold contributions (Shg/hq) and {Ahg/hg) can also 
be found in [22! • 



(A.7) 



B Higgs— Boson Couplings to Third Generation Sfermions 



Here we present the Higgs-sfermion-sfermion couplings in the weak-interaction basis. The 
couplings r"-^*-^ are given in the {f^, fn) basis by 



■p(/>l£*t 

Y4>2t*t 

■p(3r*T 
Y<t>lf*f 



_]_( ihlispAl+Cfsfx) 

y/2\-i hisf^Ab + C/3/i*) 



_]_( ih:{cpA: + sf3i^) 

V2 V -iht{cf3At + spfi*) 



— h* A* 
2„,^ I l^/2„ 



-\ht\'^vsp 



L/i* 4* 



\h,\'^vcp + \{g^-g'^)vcp -^,KK 

T^KAr -\hr\ ^VCf3 + lg%Cp 







V2 



V2' 
1 /2 

^g vsp 



[9^ + 9'^ j VC/3, 



1 

4l«' 



The coupling r^^'**'^ is given in the LR basis by 

^H+u*d _ I + \hd? - 9^)vspCfi h*a{spA*a + Cfsfi) 



K {cpA^ + spii* 
23 



g'nvsp. (B.i) 



r^^'*'^ = ^{\K\''-g')vspcp, r^^^^^- = K{s^Al + cp^i) (b.2) 
C Higgs— Boson Self— Couplings 

In the following, we list all the effective trilinear and quartic Higgs-boson self-couplings of 
the Higgs weak eigenstates, which can be expressed in terms of the conventional quartic 
couplings Ai,2,...,7 of the Higgs potential [12] obtained in an expansion of the effective Higgs 
potential up to operators of dimension 4. The trilinear couplings of the neutral Higgs 
bosons fHl are given by 



3 
3 

^010202 = C/3 A34 + cp afJeAs + 2^/3 ^eAy , 

f 0102. = -25>mA5 - spcp '^m (Ag + A7) , 

1 
2 



9^,^^a = -S/3C/3 ^^mAs - -(l + 24)5>mA7, 
9^,aa = 4^/3^1 + 4^34 - c/3(l + 4) J^eAg + ^ Sp{s} - 2cl) ^eXe + is/jc^a^eAy, 
S'^jaa = ■S/34A2 + 4-^34 - + 4) JJeAs + ^sJc/^^f^eAe + ^ c^(4 - 2^) 3^eA7 , 
9aaa = s pCp QjiiX^ " ^ ^ ^"^^6 " ^ 4 '^mXr , (C.l) 
with A34 = I (A3 + A4). The effective couplings [IE] read: 

9^^H+jj- = 2sJc/3Ai + cjA3 - sJc^A4 - 23^0^^6X5 + s^(4 - 2cJ) JReAe 
+ S/jc^JJeAy , 

3^^H+H- = 2S/34A2 + sjA3 - S/34A4 - 2s^cj3fJeA5 + s^cp^eXQ 

+ cp{cl-2sl)^eXj, 
9aH+H- = 2s/3C/3 55mA5 - s^'^mXe - S^mAy . (C.2) 

The quartic couplings for the neutral Higgs bosons [22] sue 

94,^4,^4,^4,^ ~ 4 '^-'^ ' 94^4^4^4,2 ~ 2 "''^^'^6 ' 94-^4-^4242 ~ 2 "^34 ~^ 2 "'^^'^S ! 
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^ 4>l4>2't>2't'2 2 ' 3 ,j,2cf>24>2't>2 ^ "^^ ' 

1 1 

^^i02aa = -2s(3C(3^eX5 + -sl^eXe + -cl^eX^, 

111 1 

9^,aaa = s QiiiX^ - ^ sj^/? ^mAg " ^ ^ S>mA7 , 

1 . 1 

9^2aaa = s^c/s QmX^ " - s'} '^irXq - -spcl^mXr, 

9aaaa ^ ^ H~ H+ ' 

(C.3) 

with the quartic couphng of the charged Higgs bosons being given by 

9h+h-h+h- = 4^^ + 4^2 + sJ4(A3 + A4) + 2sJ43fJeA5 - 2sJc/33?eA6 

- 2s/343fJeA7 . (C.4) 

Finally, the remaining quartic couplings involving the charged Higgs boson pairs, gapu+n^^ 
are given by 

9,,,,H+H- = + ^4^3 - w3f?eA6, 

9^-,4.2H+H- = -■5/3C/3A4 - 2spCp'^eX5 + sfj'ifieXQ + c^^eXr , 



9,,,,H+H- = + ^^As - spc^^eX 



7 ; 



9^^^^+^- = 2spCf^'^mX5 - Sf^cp'^mXe - c^'^mXj 
9^,aH+H- = 24c/3S>mA5 - s^^mXe - s^cJSJmA 



7 ■ 

9„„ n+M- 9 TT+ TT- TT+ TT- ' i.^'^^ 



D Common Blocks 

Here we list three common blocks for the SM and SUSY parameters, which are filled from 
two input arrays SMPARA_H and SSPARA_H. 

• /HC_SMPARA/: This common block is for the SM parameters. 
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COMMON /HC_SMPARA/ AEM_H , ASMZ_H , MZ_H , SW_H , ME_H , MMU_H , MTAU_H , MDMT_H 

, MSMT_H , MBMT_H , MUMT_H , MCMT_H , MTPOLE_H , GAMW_H 
, GAMZ_H , EEM_H , ASMT_H , CW_H , TW_H , MW_H , GW_H , GP_H 
,V_H,GF_H,MTMT_H 



AEM_H = acm(Mz) 
SW_H = sin Ow 
MTAU_H = rrir [GeV] 
MBMT_H = mb{m\°^'') [GeV] 
MTPOLE_H = mr'° [GeV] 
EEM_H = e = (47raem(M^))i/2 
TWJI = ianOw 
GP_H = g' = e/ cosOw 
MTMT_H = mt(mr''') 



ASMZ_H = as{Mz) 
ME_H = me [GeV] 
MDMT_H = md{m\°^'') [GeV] 
MUMT_H = m„(mr''') [GeV] 
GAMW_H = Vw [GeV] 
ASMT_H = a,(mr'^) 
MW_H = Mw = Mz cos Ow 
V_H = 2Mw/g 



MZ_H = Mz [GeV] 
MMU_H = [GeV] 
MSMT_H = m,(mr^') [GeV] 
MCMT_H = mc(mr^') [GeV] 
GAMZ_H = Tz [GeV] 
CW_H = cos 9w 
GW_H = g = e/ smOw 
GF_H = Gf 



I polc^ as{Mz) 2 

aa\mK I = jYT--, with p„, = 11 rit , 

l + /5.,^ln(mr^VM|) "^"^ 3 ^' 



2 



^pole 

mt(mr'°) = , * , pole, ■ (D.l) 

1 I (m^ ) 

/HC_RSUSYPARA/: This is for the real SUSY parameters. 

COMMON /HC_RSUSYPARA/ TB_H , CB_H , SB_H , MQ3_H , MU3_H , MD3_H , ML3_H , ME3_H 



TB_H = tan /? CB_H = cos /? SB_H = sin /? MQ3_H = Mq^ 
MU3_H = MD3_H = Mf^^ ML3_H = M^^ ME3_H = M^^ 

• /HC_CSUSYPARA/: This is for the complex SUSY parameters. 

COMPLEX* 16 MU_H , M1_H , M2_H , M3_H , AT_H , AB_H , ATAU_H 

COMMON /HC_CSUSYPARA/ MU_H,M1_H,M2_H,M3_H, AT_H, AB_H, ATAU_H 

MU_H = II [GeV] M1_H = Mi [GeV] M2_H = Ms [GeV] M3_H = M3 [GeV] 
AT_H = At [GeV] AB_H = [GeV] ATAU_H = Ar [GeV] 



E Sample Outputs 

Here we show the results of a test run of the code CPsuperH for the CPX scenario of MSSM 
Higgs-sector CP violation. 



26 



• IFLAG_H(1) = 1: The list of the SM and SUSY input parameters 



Standard Model Parameters in /HC_SMPARA/ 



AEM_H 







7812E-02 


alpha_em(MZ) 






ASMZ_H 


— 





1172E+00 


alpha_s(MZ) 






MZ_H 


— 





9119E+02 


Z boson mass in GeV 






SW_H 


_ 





4808E+00 


sinTheta_W 






ME_H 


_ 





5000E-03 


electron mass in GeV 






MMU_H 


— 





1065E+00 


muon mass in GeV 






MTAU_H 


— 





1777E+01 


tau mass in GeV 






MDMT_H 


— 





6000E-02 


d-quark mass at M_t~pole 


in 


GeV 


MSMT_H 


— 





1150E+00 


s-quark mass at M_t~pole 

T. ST 


in 


GeV 


MBMT_H 


— 





3000E+01 


b-quark mass at M_t~pole 

T. ST 


in 


GeV 


MUMT_H 


— 





3000E-02 


u-quark mass at M_t~pole 


in 


GeV 


MCMT_H 


- 





6200E+00 


c-quark mass at M_t~pole 


in 


GeV 


MTPOLE_H 


— 





1750E+03 


t-quark pole mass in GeV 






GAMW_H 


— 





2118E+01 


Gam_W in GeV 






GAMZ_H 


= 





2495E+01 


Gam_Z in GeV 






EEM_H 







3133E+00 


e = (4*pi*alpha_em)"l/2 






ASMT_H 







1072E+00 


alpha_s(M_t"pole) 






CW_H 







8768E+00 


cosTheta_W 






TW_H 







5483E+00 


tanTheta_W 






MW_H 







7996E+02 


W boson mass MW = MZ*CW 






GW_H 







6517E+00 


SU(2) gauge coupling gw= 


e/ s 


.W 


GP_H 







3573E+00 


U(1)_Y gauge coupling gp= 


=e/c 


:_W 


V_H 







2454E+03 


V = 2 MW / gw 






GF_H 







1174E-04 


GF=sqrt(2)*gw~2/8 MW"2 in GeV"2 


MTMT_H 







1674E+03 


t-quark mass at M_t"pole 


in 


GeV 


Real SUSY 


Parameters in /HC_RSUSYPARA/ 






TB_H 







5000E+01 


tan(beta) 






CB_H 







1961E+00 


cos (beta) 






SB_H 







9806E+00 


sin(beta) 






MQ3_H 







5000E+03 


M_tildeQ_3 in GeV 






MU3_H 







5000E+03 


M_tildeU_3 in GeV 






MD3_H 







5000E+03 


M_tildeD_3 in GeV 






ML3_H 







5000E+03 


M_tildeL_3 in GeV 






ME3_H 







5000E+03 


M_tildeE_3 in GeV 
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Complex SUSY Parameters in /HC_CSUSYPARA/ 



1 MU_H 1 




0. 


.2000E+04: 


:Mae;. 


of 


MU 


parameter 


in 


GeV 


|M1_H| 


_ 


0. 


.5000E+02: 


:Maff. 


of 


Ml 


parameter 


in 


GeV 


1 M2_H 1 


_ 


0. 


.lOOOE+03: 


:Maff. 

O 


of 


M2 


parameter 


in 


GeV 


1 M3_H 1 


_ 


0. 


. lOOOE+04: 


:Maff. 

O 


of 


M3 


parameter 


in 


GeV 


1 AT_H 1 


= 


0. 


. lOOOE+04: 


:Mag. 


of 


AT 


parameter 


in 


GeV 


1 AB_H 1 




0. 


. lOOOE+04: 


:Mag. 


of 


AB 


parameter 


in 


GeV 


1 ATAU_H 1 




0. 


. lOOOE+04: 


:Mag. 


of 


ATAU parameter : 


Ln GeV 


ARG(MU_H) 




0. 


.OOOOE+00: 


:Arg. 


of 


MU 


parameter 


in 


Degree 


ARG(M1_H) 




0. 


.OOOOE+00: 


:Arg. 


of 


Ml 


parameter 


in 


Degree 


ARG(M2_H) 




0. 


.OOOOE+00: 


:Arg. 


of 


M2 


parameter 


in 


Degree 


ARG(M3_H) 




0. 


,9000E+02: 


:Arg. 


of 


M3 


parameter 


in 


Degree 


ARG(AT_H) 




0. 


.9000E+02: 


:Arg. 


of 


AT 


parameter 


in 


Degree 


ARG(AB_H) 




0. 


.9000E+02: 


:Arg. 


of 


AB 


parameter 


in 


Degree 


ARG(ATAU_H) 




0. 


.9000E+02: 


:Arg. 


of 


ATAU parameter : 


Ln Degree 



Charged Higgs boson pole mass : 0.3000E+03 GeV 



• 1FLAG_H(2) = 1: The masses and mixing matrix of the neutral Higgs bosons 



Masses and Mixing Matrix of Higgs bosons : 

HMASS_H(1) and OMIX_H(A,I) 



HI Pole Mass = 0.1188E+03 GeV 

H2 Pole Mass = 0.2703E+03 GeV 

H3 Pole Mass = 0.2981E+03 GeV 

Charged Higgs Pole Mass = 0.3000E+03 GeV [SSPARA_H(2)] 

[HI] [H2] [H3] 

[phi_l] / 0.2451E+00 -.3373E+00 -.9089E+00 \ 
0(IA,IH)=[phi_2] I 0.9694E+00 0.7532E-01 0.2335E+00 I 
[ a ] \ -.1030E-01 -.9384E+00 0.3454E+00 / 



• IFLAG_H(4) = 1: The masses and mixing matrices of the charginos and neutrahnos 



Chargino Masses and Mixing Matrices : 
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MC_H(I), UL_H(I,A), and UR_H(I,A) 


MCI = 0.9861E+02 GeV 


MC2 = 0.2003E+04 GeV 


UL_H = 




/(0.9984E+00 O.OOOOE+00) 


(-.5603E-01 O.OOOOE+00) \ 


\(0.5603E-01 O.OOOOE+00) 


(0.9984E+00 O.OOOOE+00) / 


UR_H = 




/(0.9999E+00 O.OOOOE+00) 


(-.1385E-01 0.5460E-09) \ 


\(0.1385E-01 O.OOOOE+00) 


(0.9999E+00 O.OOOOE+00) / 



Neutralino Masses MN_H(I) and Mixing Matrix N_H(I,A) 



MNl = 


4960E+02 GeV 




MN2 = . 9862E+02 


GeV 


MN3 = 


2001E+04 GeV 




MN4 = . 2003E+04 


GeV 


N_H(1,1) 


= (0 


9996E+00 





OOOOE+00) 




N_H(1,2) 




1462E-01 





OOOOE+00) 




N_H(1,3) 


= (0 


2218E-01 





OOOOE+00) 




N_H(1,4) 




4962E-02 





OOOOE+00) 




N_H(2,1) 




1553E-01 





OOOOE+00) 




N_H(2,2) 




9991E+00 





OOOOE+00) 




N_H(2,3) 


= (0 


3931E-01 





OOOOE+00) 




N_H(2,4) 




9704E-02 





OOOOE+00) 




N_H(3,1) 


= (0 


OOOOE+00 




1186E-01) 




N_H(3,2) 


= (0 


OOOOE+00 





2112E-01) 




N_H(3,3) 


= (0 


OOOOE+00 





7066E+00) 




N_H(3,4) 


= (0 


OOOOE+00 





7072E+00) 




N_H(4,1) 


= (0 


1867E-01 





OOOOE+00) 




N_H(4,2) 




3494E-01 





OOOOE+00) 




N_H(4,3) 




7062E+00 





OOOOE+00) 




N_H(4,4) 


= (0 


7069E+00 





OOOOE+00) 




IFLAG_H(5) = 1: 


The couplin; 


2;s of the hghtest Higg 


s boson 



The Lightest Higgs H_l Couplings : NHC_H(NC,1) 



HI 


e e 


[NC= 1] 


GF= 


cO 


2038E-05 





OOOOE+00) 






[NC= 2] 


GS= 


^0 


1250E+01 





OOOOE+00) 






[NC= 3] 


GP= 


^0 


5148E-01 





OOOOE+00; 


HI 


mu mu 


[NC= 4] 


GF= 


[0 


4340E-03 





OOOOE+00) 






[NC= 5] 


GS= 




1250E+01 





OOOOE+00) 






[NC= 6] 


GP= 




1 — A A "I — ' A 

5148E-01 





OOOOE+00) 


HI 


tau tau 


[NC= 7] 


GF= 


Co 


7242E-02 





OOOOE+00) 






[NC= 8] 


GS= 


[0 


1250E+01 





OOOOE+00) 






[NC= 9] 


GP= 


[0 


5148E-01 





OOOOE+00) 


HI 


d d 


[NC=10] 


GF= 


[0 


2445E-04 





OOOOE+00) 






[NC=11] 


GS= 


[0 


1250E+01 





OOOOE+00) 






LNC=12J 


GP= 


io 


5148E-01 





OOOOE+00) 


HI 


s s 


[NC=13] 


GF= 


[0 


4687E-03 





OOOOE+00) 






[NC=14] 


GS= 


'0 


1250E+01 





OOOOE+00) 






[NC=15] 


GP= 


CO 


5148E-01 





OOOOE+00) 


HI 


b b 


[NC=16] 


GF= 


Co 


1223E-01 





OOOOE+00) 






[NC=17] 


GS= 


CO 


1246E+01 





OOOOE+00) 






LNC=18J 


GP= 




1741E-01 





OOOOE+00) 


HI 


u u 


[NC=19] 


GF= 


CO 


1223E-04 





OOOOE+00) 






[NC=20] 


GS= 


Co 


9886E+00 





OOOOE+00) 






LNC=21J 


GP= 


Co 


2059E-02 





OOOOE+00) 


HI 


c c 


LNC=22J 


GF= 


Co 


2527E-02 





OOOOE+00) 






[NC=23] 


GS= 


CO 


9886E+00 





OOOOE+00) 






[NC=24] 


GP= 


CO 


2059E-02 





OOOOE+00) 


HI 


t t 


[NC=25] 


GF= 


CO 


6821E+00 





OOOOE+00) 






[NC=26] 


GS= 


CO 


9892E+00 





OOOOE+00) 






[NC=27] 


GP= 


CO 


4501E-02 





OOOOE+00) 


HI 


Nl Nl 


LNC=28J 


GF= 


CO 


3258E+00 





OOOOE+00) 






[NC=29] 


GS= 




5767E-02 





OOOOE+00) 






[NC=30] 


GP= 


CO 


1317E-03 





OOOOE+00) 


HI 


N2 N2 


[NC=31] 


GF= 


CO 


3258E+00 





OOOOE+00) 






[NC=32] 


GS= 




1886E-01 





OOOOE+00) 






[NC=33] 


GP= 


CO 


4125E-03 





OOOOE+00) 


HI 


N3 N3 


[NC=34] 


GF= 


CO 


3258E+00 





OOOOE+00) 






[NC=35] 


GS= 


CO 


1415E-01 





OOOOE+00) 






[NC=36] 


GP= 


CO 


1576E-03 





OOOOE+00) 


HI 


N4 N4 


[NC=37] 


GF= 


CO 


3258E+00 





OOOOE+00) 







[NC=38] 


GS= 


:o 


3878E-01 





OOOOE+00) 






[NC=39] 


GP= 




3866E-03 





OOOOE+00) 


HI 


Nl N2 


[NC=40] 


GF= 


::o 


3258E+00 





OOOOE+00) 






[NC=41] 


GS= 




1043E-01 





OOOOE+OO) 






[NC=42] 


GP= 


:o 


2331E-03 





OOOOE+OO) 


HI 


Nl N3 


[NC=43] 


GF= 


[0 


3258E+00 





OOOOE+00) 






[NC=44] 


GS= 




1602E-02 





OOOOE+00) 






[NC=45] 


GP= 


^0 


1443E+00 





OOOOE+00) 


HI 


Nl N4 


[NC=46] 


GF= 


::o 


3258E+00 





OOOOE+00) 






[NC=47] 


GS= 




2413E+00 





OOOOE+00) 






[NC=48] 


GP= 




2402E-02 





OOOOE+00) 


HI 


N2 N3 


[NC=49] 


GF= 




3258E+00 





OOOOE+00) 






[NC=50] 


GS= 




2820E-02 





OOOOE+00) 






[NC=51] 


GP= 


:o 


2540E+00 





OOOOE+OO) 


HI 


N2 N4 


[NC=52] 


GF= 




3258E+00 





OOOOE+00) 






[NC=53] 


GS= 


[0 


4247E+00 





OOOOE+00) 






[NC=54] 


GP= 




4228E-02 





OOOOE+00) 


HI 


N3 N4 


[NC=55] 


GF= 




3258E+00 





OOOOE+00) 






[NC=56] 


GS= 




2470E-03 





OOOOE+00) 






[NC=57] 


GP= 




2789E-03 





OOOOE+OO) 


HI 


C1+ ci- 


[NC=58] 


GF= 


:o 


4608E+00 





OOOOE+00) 






[NC=59] 


GS= 




2714E-01 





OOOOE+OO) 






[NC=60] 


GP= 


'0 


5936E-03 




2135E-18) 


HI 


C1+ C2- 


[NC=61] 


GF= 


'0 


4608E+00 





OOOOE+OO) 






[NC=62] 


GS= 


^0 


6058E+00 





4035E-02) 






[NC=63] 


GP= 




6042E-02 




3618E+00) 


HI 


C2+ Cl- 


[NC=64] 


GF= 


::o 


4608E+00 





OOOOE+00) 






[NC=65] 


GS= 


'0 


6058E+00 




4035E-02) 






[NC=66] 


GP= 




6042E-02 





3618E+00) 


HI 


C2+ C2- 


[NC=67] 


GF= 




4608E+00 





OOOOE+00) 






[NC=68] 


GS= 


:o 


5771E-01 





OOOOE+00) 






[NC=69] 


GP= 




2527E-03 





7454E-19) 


HI 


V V 


[NC=70] 


G = 


CO 


9987E+00 





OOOOE+00) 


HI 


STl* STl 


[NC=71] 


G = 


[0 


2184E+01 




3755E-16) 


HI 


STl* ST2 


[NC=72] 


G = 


[0 


2447E+00 




1144E+00) 


HI 


ST2* STl 


[NC=73] 


G = 


[0 


2447E+00 





1144E+00) 


HI 


ST2* ST2 


[NC=74] 


G = 




3987E+01 




4869E-17) 


HI 


SBl* SBl 


[NC=75] 


G = 


[0 


4330E+00 





1506E-18) 


HI 


SBl* SB2 


[NC=76] 


G = 


CO 


2070E-01 





8020E-02) 


HI 


SB2* SBl 


[NC=77] 


G = 


CO 


2070E-01 




8020E-02) 



HI 


SB2* SB2 


[NC=78] 


G 






5584E+00,0 


1683E-18) 


HI 


STAl* STAl 


[NC=79] 


G 





^0 


2300E+00,- 


1144E-17) 


HI 


STAl* STA2 


[NC=80] 


G 


_ 


:o 


1602E-02,0 


6836E-02) 


HI 


STA2* STAl 


[NC=81] 


G 


_ 


:o 


1602E-02,- 


6836E-02) 


HI 


STA2* STA2 


[NC=82] 


G 


— 




3549E+00,- 


1785E-17) 


HI 


SNU3* SNU3 


[NC=83] 


G 


_ 


:o 


1246E+00,0 


OOOOE+00) 


HI 


glue fflue 


[NC=84] 


S 





^0 


5827E+00,0 


3665E-01) 






[NC=85] 


P 


_ 




5195E-02,- 


5135E-03) 


HI 


CH+ CH- 


[NC=86] 


G 


— 




4047E-01,0 


OOOOE+00) 


HI 


CH+ W- 


[NC=87] 


G 


= 


:;- 


5024E-01,0 


1030E-01) 


HI 


photon photon 


[NC=88] 


S 






6557E+01,0 


2443E-01) 






[NC=89] 


P 






1385E-01,- 


3423E-03) 


HI 


glue glue 


(M=0) [NC=90] 


S 






1427E+01,0 


1915E-17) 






[NC=91] 


P 






1291E-01,0 


OOOOE+00) 


HI 


photon photon 


(M=0) [NC=92] 


S 






4878E+01,0 


5235E-17) 






[NC=93] 


P 






1728E-02,- 


4811E-18) 



• IFLAG_H(6) = 1: The decay width and branching fractions of the Ughtest Higgs boson 



Neutral Higgs Boson Decays with ISMN = 50 : ISUSYN = 50 



DECAY MODE [ IM] WIDTH [GeV] BR[SM] BR [TOTAL] 



HI 


-> 


e 


e [ 


1] : 





3070E-10 





5766E-08 





5760E-08 


HI 


-> 


mu 


mu [ 


2] : 





1393E-05 





2616E-03 





2613E-03 


HI 


-> 


tau 


tau [ 


3] : 





3873E-03 





7274E-01 





7266E-01 


HI 


-> 


d 


d [ 


4] : 





1686E-07 





3166E-05 





3163E-05 


HI 


-> 


s 


s [ 


5] : 





6193E-05 





1163E-02 





1162E-02 


HI 


-> 


b 


b [ 


6] : 





4163E-02 





7820E+00 





7812E+00 


HI 


-> 


u 


u [ 


7] : 





2632E-08 





4944E-06 





4939E-06 


HI 


-> 


c 


c [ 


8] : 





1124E-03 





2111E-01 





2109E-01 


HI 


-> 


t 


t [ 


9] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


W 


W [ 


10] : 





4106E-03 





7711E-01 





7703E-01 


HI 


-> 


Z 


Z [ 


11] : 





3303E-04 





6203E-02 





6197E-02 


HI 


-> 


HI 


Z [ 


12] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


H2 


Z [ 


13] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


HI 


HI [ 


14] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


HI 


H2 [ 


15] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 
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HI 


-> 


H2 


H2 


: 16] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


ph 


ph 


: 17] : 





9427E-05 





1771E-02 





1769E-02 


HI 


-> 


gl 


gl 


: 18] : 





2004E-03 





3763E-01 





3760E-01 


HI 


TOTAL (SM) 


: 50] : 





5324E-02 





lOOOE+01 





9990E+00 


HI 


-> 


Nl 


Nl 


: 51] : 





5557E-05 





1044E-02 





1043E-02 


HI 


-> 


Nl 


N2 


: 52] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


Nl 


N3 


: 53] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


Nl 


N4 


: 54] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


N2 


N2 


: 55] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


N2 


N3 


: 56] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


N2 


N4 


: 57] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


N3 


N3 


: 58] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


N3 


N4 


: 59] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


N4 


N4 


: 60] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


C1+ 


ci- 


: 61] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


C1+ 


C2- 


: 62] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


C2+ 


ci- 


: 63] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


C2+ 


C2- 


: 64] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


STl* 


STl 


: 65] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


STl* 


ST2 


: 66] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


ST2* 


STl 


: 67] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


ST2* 


ST2 


: 68] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


SBl* 


SBl 


: 69] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


SBl* 


SB2 


: 70] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


SB2* 


SBl 


: 71] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


-> 


SB2* 


SB2 


: 72] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


->STA1*STA1 


731 ■ 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


->STA1*STA2 


: 74] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


->STA2*STA1 


: 75] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


->STA2*STA2 


: 76] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


->SNU3*SNU3 


: 77] : 





OOOOE+00 





OOOOE+00 





OOOOE+00 


HI 


TOTAL (SUSY) 


:i00] : 





5557E-05 





1044E-02 





1043E-02 


HI 


TOTAL 




:ioi] : 





5330E-02 





lOOlE+01 





lOOOE+01 



* Note : WIDTH=GAMBRN(IM,1,1) , BR[SM] =GAMBRN(IM,2, 1) 

and BR [TOTAL] =GAMBRN(IM, 3,1) 
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Table 4: The couplings of the lH-th neutral Higgs boson to two particles specified with the 
index NC, NHC_H(NC, IH) . For the definitions of the couplings to two fermions, gj, g^ , and 
g^, see Eq. ^ and TableU\ 



NC 


Coupling 


NC 


Coupling 




NC 


Coupling 


y 


NC 


Coupling 


1 




26 






51 


^£hx!^xU' 


^k^xlxl 


76 


9Hjibfb2 


2 


s 

9Hme+e- 


27 


SUmtt 




52 


^X« 




77 


dHmbfbi 


3 




28 


9x° 




53 


s 

^H^^xlxV 


s 

^H^^xlxl 


78 


9Hjibfb2 


4 




29 






54 


P 

^H^KXlxV 


P 

3h,^xIx°2 


79 


9Hmr(fi 


5 


s 

9Hl^^l+^l- 


30 


p 




55 


^X« 




80 


9Hmf*f2 


6 


p 


31 


^XO 




56 


s 

^H^^xlxV 


s 

3h^^xIxI 


81 


9Hmf*Ti 


( 




OZ 


^^IHX§X^ 




b ( 


^H.^x'ixl' 


^H^^xlxl 


OZ 




P 

o 


yHmr+T- 


OO 






OO 


9x± 




OO 


9HiiiU*I>t 


Q 




0*1 


9x° 






^^^iHX+xr 








10 


gd 


35 


s 

^^^iHxgxi 




60 


p 

fl'rr -+-- 
^HmXiXi 




85 




11 


s 

dHmdd 


36 
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61 


9x^ 




86 




12 


p 

duiudd 


37 


^XO 




62 
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9ij -+ — 




87 




13 


9s 


38 


^^i„x2x2 




63 


^Hmxtx2 




88 


SUMhJ 


14 




39 
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89 


pUMhJ 


15 




40 
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9b 


41 
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9 Hjj^VV 
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dHjYiUU 
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Table 5: The trilinear (\iC=l-13) and quartic (iiC=14-35) Higgs-boson self-couplings, 
SHC_H(NC). We note that SHC_H(10+IH)=NHC_H(86 , IH) for m=l-3. 



NC 


Coupling 


NC 


Coupling 


NC 


Coupling 


NC 


Coupling 


1 




11 


9 H^H+H- 


21 


9H3H2H2H1 


31 


9 H-jHiH+H- 
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9H3H3H2 


12 


9 H2H+H- 
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9H-jH2HiHi 


32 


9h2 H2H+H- 
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9h.jH.jHj 
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9 H-^H+H- 
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9 H^HiHiHi 
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9 H2 H2 H2H2 
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15 


9h^H^H^H2 
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Table 6: The charged Higgs-boson couplings to two particles specified with the index NC, 
CHC_H(NC). For the definitions of the couplings to two fermions, Qff, , and g^, see 
Eq. ^ and Table\^ 
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Table 7: The decay mode with the index IM of the neutral Higgs bosons Hm used in 
GAMBRN(IM,IWB,IH). When IWB=i, IM=ISMN+ISUSYN+1 is for the total decay width of the 
neutral Higgs boson Hm and IM=ISMN and IM=ISMN+ISUSYN for the subtotal decay widths 
into SM particles fGAMSM^ and into SUSY particles (GMSVSY), respectively. Note that the 
current version o/ CPsuperH does not compute (loop-induced) absorptive phases, i.e., it 
currently returns equal decay widths into CP-conjugate final states. 



IM 



Decay Mode 



IM 


Decay Mode 


11 


Hm - 


ZZ 


12 


Hm - 


HiZ 


13 


Hm - 


-^H^Z 


14 


Hm - 


--H,H, 


15 


Hm - 


H1H2 


16 


Hm - 


H2H2 


17 


Hm - 




18 


Hm - 


^99 



IM 



Decay Mode 



1 


Hm 


ee 


2 


Hm — > 


fifi 


3 




rf 


4 


Hm — > 


dd 


5 


Hm — * 


ss 


6 


Hm 


hi 


7 


Hm —> 


uu 


8 


Hm —> 


cc 


9 


Hm 


tt 


10 


Hm 


WW 



ISMN 



GAMSM 



IM 



Decay Mode 



IM 


Decay Mode 


ISMN+11 


Hm - 


xtxi 


ISMN+12 


Hm - 


^ X1X2 


ISMN+13 


Hm - 


X2X1 


ISMN+14 


Hm - 


X2X2 


ISMN+15 


Hm - 


^t\ti 


ISMN+16 


Hm - 




ISMN+17 


Hm - 


t^ti 


ISMN+18 


Hm ~ 


H> ^2^2 


ISMN+19 


Hm - 


-.bib. 


ISMN+20 


Hm - 


-.blb2 



IM 



Decay Mode 



ISMN+1 
ISMN+2 
ISMN+3 
ISMN+4 
ISMN+5 
ISMN+6 
ISMN+7 
ISMN+8 
ISMN+9 
ISMN+10 



Hm 
Hm 
Hm 
Hm 
Hm 
Hm 
Hm 
Hm 
Hm 
Hm 



xlxl 

x'lXl 

xlxl 

X'2X'2 
X'2Xl 

xtxl 
xlxl 
xtxl 
xlxl 



ISMN+21 
ISMN+22 
ISMN+23 
ISMN+24 
ISMN+25 
ISMN+26 
ISMN+27 



Hm 
Hm 
Hm 
Hm 
Hm 



b;b2 
T1T2 

T*Ti 
T*f2 



ISMN+ISUSYN GAMSUSY 
ISMN+ISUSYN+1 GAMSM+GAMSUSY 



41 



Table 8: The decay mode with the index IM of the charged Higgs boson used in 
GAMBRC(IM,IWB). When IWB=i, IM=ISMC+ISUSYC+1 is for the total decay width of the 
charged Higgs boson and IM=ISMC and IM=ISMC+ISUSYC for the subtotal decay widths into 
SM particles (GMSH) and into SUSY particles fGAMSUSY^, respectively. 
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Figure 1: The branching ratios and total decay widths of the MSSM Higgs bosons, taking into 
account only the decays into SM particles. All the CP phases are set to zero and tan /5 = 1.5 
is taken. For the comparison with HDECAY, the threshold corrections are not included and 
we assume the 'maximal mixing' scenario: lA.fe.rl = V^A^susy with the common SUSY 
scale MsusY = Mq^ = M^^ = M^^ = = M^^ = 1 TeV and \fi\ = 100 GeV. The results 
are consistent with those from the code HDECAY. 



43 



BR(H,) 



10 



-1 



10 



10 



-3 





■ bb 




^ XT 




t 




g§// WW 




I ..YY 

1 1 1 1 / 1 1 



60 80 

BR(H^) 



100 120 

MjjJGeV] 



10 



-1 



10 



10 



HI HI 


WW mm 

X 1 * 

f 


tt^ 

/ V 
/ ^ 




. XT' 


^ - . 

ZZ 


\ 

\ 




f g /' 
-'0'' 

L \ 




V ^ 


: \ 

' i 


l\ * 
l\ * 

1 \ ^ 

\ ^ 
1 I 


t 

* 1 ^ 


1 " "1 



10 



5x10 
Mjj3[GeV] 



BR(HJ 



10 



-1 



10 



-3 



10 



10 



bb 



XT 



10 

r(H) [GeV] 



10 



10 



10 



-1 



10 



HI Z 



gg/ 



5x10 



Mjj2[GeV] 













H3 


/ t 


- * 
1 




H2 \.'J 




^ 1 




1 1 


1 

f 

1 


1 1 1 



5x10 
M^LGeV] 



Figure 2: The branching ratios and total decay widths of the MSSM Higgs bosons, taking into 
account only the decays into SM particles. All the CP phases are set to zero and tan (3 = 4 is 
taken. We assume the CPX scenario: \M-i\ = 1 TeV, |/i| = 4Msusy O'lT'd |At,b,r| = 2Msusy 
with the common SUSY scale MgusY = Mr, = Mfj = = = = 0.5 TeV. 
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Figure 3: The same as in Fig. but with non-trivial CP phases: = = 
$3 = 90°. 
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Figure 4: The dependence of the branching ratios of Higgs bosons into superparticles on the 
phase 0^, for tan/3 = 5, Mf,f = 0.3 TeV, M^^ = M^^ = M^^ = M^.^ = M^^ = 0.5 TeV, 
= 250 GeV, \Mi\ = 50 GeV, IM2I = 150 GeV, {M^l = 0.5 TeV, \At\ = \Ab\ = \Ar\ = 1.2 
TeV, ^2 = ^3 = 0, and = = = The left frame shows results for Hi 

for several choices of ^i; in this case the only contributing final state is XiXi- The right 
frame shows results for the heavier Higgs bosons, where the solid (dotted) lines are for 
$1 = (180°); in this case heavier neutralinos and charginos contribute, but no sfermions. 
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